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Abstract. - We discuss with the aid of random walk arguments and exact numerical com- 
putations the magnetization properties of one-dimensional random field chains. The ground 
state structure is explained in terms of absorbing and non-absorbing random walk excursions. 
At low temperatures, the magnetization profiles follow those of the ground states except at 
regions where a local random field fluctuation makes thermal excitations feasible. This follows 
also from the non-absorbing random walks, and implies that the magnetization length scale is 
a product of these two scales. It is not simply given by the Imry-Ma-like ground state domain 
size nor by the scale of the thermal excitations. 



In statistical mechanics of random systems the search for universality can be interpreted 
geometrically. That is, if the introduction of disorder into a system is relevant, the real-space 
properties of the physical states can be understood in terms of scaling exponents. These 
describe the fluctuations of a domain wall, or the behavior of a spin-spin correlation function. 
The central ingredient is that the configurational energy is coupled to geometric fluctuations. 
Consider a domain wall in a magnet. If the spatial fluctuations are described by a roughness 
exponent £, then there is an associated exponent describing the free or ground state energy 
fluctuations. Assuming that the 'zero temperature fixed point' scenario is true or that the 
entropy is irrelevant at low enough temperatures, this is all what is needed to describe the 
physics. The system evolves via Arrhenius-like dynamics so that the cost of moving in the 
energy landscape is given by the usual exponential factor exp(A_E/3), where (3 — 1/T and T 
is the temperature, and AE ~ ft relates the cost to the scale length of the perturbation I. 

Consider now a random magnet. It has a ground state (GS) which is described exactly 
by the positions and arrangement of the domain walls. Examples abound in particular in 
Ising systems, where non-trivial GSs exist for spin glasses and random field systems jjj. In 
this work we investigate with random walk arguments and exact numerical computations 
how the aforementioned picture applies in the case of one-dimensional random field chains. 
We find that for arbitrary field distributions M the GS structure can be understood via the 
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Fig. 1 - Example for a ground state configuration (bold line) of a chain segment for a particular real- 
ization of random fields as compared to the equilibrium expectation values of the local magnetization 
(<7j) at temperature T = 0.45J. In the leftmost domain nearly no melting occurs, the region right of 
i — 125 is a typical example for boundary melting and in the rightmost domain the bulk melting is 
so strong that even magnetization reversal occurs. 



random walk picture, which is compared to exact numerical GS computations. At finite 
temperatures we resort to scaling arguments based on this random walk picture, and again to 
an exact numerical determination of the magnetization. This allows us to make conclusions 
about the behavior in the same sample in both cases. Our main finding is, in addition to 
recovering the GS from the random walk picture, the emergence of two relevant length scales. 
These arise from the zero-temperature length scale of the domains and the typical size of 
'easy' excitations at a given temperature. The latter changes the correlation length of the 
magnetization, and thus leads to the fact that in our case the low-temperature physics is 
characterized not only by the zero-temperature scaling. The Id RFIM has received recently 
attention || since it is simple enough that decimation-type real space renormalization can be 
applied to domain wall dynamics (each DW undergoes logarithmic Sinai diffusion which 
can be compared with our findings concerning the asymptotic state of such processes. The 
zeroes of the magnetization profile simply denote the equilibrium positions of domain walls at 
T > 0, and the extra physics consists of additional domain walls added to the GS structure. 
The chain is also the simplest example of a random magnet with a competition of non-trivial 
GS and thermal excitations (e.g. random bond Ising magnets have a trivial GS). 
In the following we investigate the Hamiltonian 



N 



H = - } j Jo- 4 cr i+ i + 2_j hi(Ji 



(1) 



i=l 



where the cr's are spins located at sites i of the chain, and hi are random fields picked from a 
suitable probability distribution P{hi) with zero mean and variance h r . For a binary distri- 
bution hi = ±hi^ r the model is equivalent to a spin glass chain (with couplings Ji = ± J) in a 
homogeneous external field h r [|[|(|. Fig. |l| shows typical GS and finite temperature (T = 0.4) 
magnetization profiles obtained from the numerical procedures described below. The GS do- 
main size is often thought to be given by the Imry-Ma argument 0, which states that the 
domain field energy balances the cost from the domain walls on a scale [l) av ~ 1/h? in ID and 
[. . .] av denotes the disorder average. This reasoning omits the global optimization behind the 
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Fig. 2 - Terminology used for the description of the RW arguments. For further details see text. 

GS; later we discuss the exact way the optimization becomes visible in. At finite but small 
temperatures the magnetization changes due to two reasons. The GS domain walls fluctuate, 
and thus the m(:r)-profile is smoothed out around the GS positions. More interestingly, there 
are regions inside domains where the magnetization can even undergo a local reversal. Both 
of the cases arise from the local random field configurations as we now demonstrate. 

The starting point for the random walk argument is the fact that any sequence S of lattice 
sites i with | J2ies hi\>2J evidently leads to a GS spin structure with <Ji = +1 \/i s S if 
Sies hi — (and <7j = — 1 Vi G S if X))es^' — — ^J) independent of the local fields hj 
at sites j S. The system can thus be split up into such absorbing excursions and into the 
remaining lattice sites, which make up so-called non-absorbing excursions. 

Fig. |^ illustrates these concepts. An absorbing excursion is a sequence of spins starting at 
some lattice site i and ending at the lattice site j > i, with the field-sum | X)iG5 ^i\ for the 
first time becoming greater or equal to 2J: 

J k 

| hi\ > 2J and | ^ h t \ < 2J Vi<k<j. (2) 

l—i l—i 

In Fig. U the left- and rightmost sequences are absorbing excursions, of length l ae . A sequence 
S' of spins from i to j > i is a non- absorbing excursion if 

j k 

ct^Z hi - and < hi <2J ^ i < k <j ( 3 ) 

l—i l—i 

where a = ±1 is the orientation of the spins within the preceding absorbing excursion. The 
length of a non-absorbing excursion is l nae . A simple 'step down' (like from Si to £2) is 
included in this definition. 

The GS now follows as a sequence of absorbing and non-absorbing excursions. It, and 
the Zeeman energy and mean domain-length can be determined with the three rules: (1) 
determine an absorbing excursion Sq for a given field configuration. If it starts at site iq, ends 
at jo, and a is the sign of its field-sum, then =~o for all k 6 So- (2) start from jo + 1 and 
find all n nae of non-absorbing excursions until the next absorbing excursion 1S1 (from i\ to ji) 
is found, whose field-sum is by definition opposite in sign to the preceding one. The sites k 
belonging to the non-absorbing excursions have the same orientation (jfc = a as those within 
iSo- The orientation of the spins at sites I within Si is opposite to the latter one, o~i = —a. 
(3) starting again at j\ + 1 the search (2) for the next absorbing excursion then leads to 



4 




Fig. 4 - Probability distribution of the domain lengths Id- Apart from a non- exponential tail which 
might be due to finite size effects the decay is exponential with decay rate v. Inset: The decay rate 
v as a function of h r . For fc r <l the data are compatible with v ~ h% (bold line). 



the overall GS. These steps actually define a fast algorithm for finding the GS, though for 
historical reasons we have used the mapping to the max-flow/min-cut problem JlJ. 
The mean domain length [ld]a, v is given by 

[ld]av{h r ) — [ %ae]av[^nae]av H~ [^aejav- (4) 

Thus the domain length consists of two distinct contributions and we need to estimate the 
/^-dependence of [n„ Qe ] av , [i noe ] a v an d ['oe]av The fieldsum over the local fields of a non- 
absorbing excursion is a RW with absorbing boundaries at J2 na e ^« ~ and J2 nae = ^ 
and random step size with zero mean and variance h r . Rescaling the step size h[ = hi/h r — > 1 
this becomes a Id RW starting from x = 1 at t = with random step lengths with mean 
zero and variance one (for a binary distribution hi — ±/i r it yields a conventional lattice 
random walk with h! i = ±1) and absorbing boundaries at x — and x = L = 2J/h r . The 
probability Pa(t, L) to be absorbed at x = within the time interval [t, t + dt] without having 
been absorbed at x — L reads || Po(t, L) oc t~ 3//2 for 1 < t < L 2 and decays exponentially 
for t > L 2 . Integration over Po(t, L) leads to [i na e]av ~ L °t h~ . The mean number [n„ e]av 
of consecutive non-absorbing excursions follows due to the fact that the probability for an 
excursion to be absorbing grows as p ae ~ 1/L oc h r |)|. Thus P(n nae ) ~ (1 — Pae) nna ° decays 
exponentially. As a consequence n nae ~ l/ln(l — bh r ) and the mean length of an absorbing 
excursion grows like h~ 2 . Finally eq. (|J) reads 

[Mav(M - ^ ln(1 Q _ &M +£l - ^ for hr^O (5) 

where one expects b < 1, a < and —a » c. Note that a//i r ln(l — 6/i r ) ~ ft,,: 2 for h r — > 
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and for h r < 0.5 no significant difference between a/[h r ln(f — b h r )] and h~ 2 can be observed. 
The asymptotic limit follows the Imry-Ma scaling, though the physics is more complicated. 

This result is confirmed by computations of exact ground states. The data was obtained 
for a Gaussian random field distribution with zero mean and variance h r and averaged over 
10 5 disorder configurations. The system size is large enough (L = 5000) such that [^] av <C L 
even for the smallest field strength h r . Fig. || shows our numerical result for the average length 
[Id — 2] av of the GS domains as a function of the field amplitude h r . In the limit h r — > oo 
I'd] av - ► 2 since all the spins align with their local fields. In the limit h r <C J the data fit well 
to the predicted form eq. (||), scaling as h~ 2 for h r — > 0. Moreover, as can be seen in Fig. [| 
the probability distribution of the domain sizes decays exponentially, with a decay rate v that 
scales inversely proportional to [ijavj i-e. v(h r ) oc h 2 . 

The field energy of a domain can be computed as a function of h r and Id by noting that 
both a single absorbing excursion and all of the non-absorbing excursions contribute. The 
former contributes a constant (2J), depending neither on h r nor on Id- Each non-absorbing 
excursion adds an amount of 0{h r ) so that the sum self-averages. The contribution of a 
single non-absorbing excursion equals £j — £j_i ~ h r , i.e. the step width of the RW. Thus 
the field energy results from the number of non-absorbing excursions in a domain, n nae , plus 
2J. From @ we learn that in the limit h r — > oo the contribution of the absorbing and 
non- absorbing walks to [ld]av scale similarly such that we expect that for a fixed domain size 

Kae('(i)]av OC ld/[lnae}a.v OC l d h r . Thus 

[E f (l d )U = 2J+[n noe (Z d )] av /i r = 2J + dh 2 r [Z d ] av (6) 

The numerics confirms this result: Fig. || shows the data for the mean Zeeman energy 
[Ef(ld)]av of domains of length Id, From the slopes of the straight lines we learn that [Ef(ld)] av 
is linear in the domain length and from the offsets that it grows like h 2 , independent of the 
field distribution P(hi). Note that from a naive random walk picture one would expect 

1 /2 

[Ef(ld)]a, v oc l d h r , which is incorrect. 

We now turn our attention to equilibrium configurations, i.e. the local magnetization 
m(x) and the domain structure at T > 0. Using numerical transfer matrix methods [10,11] to 
compute the partition function Zjv we can compute the exact expectation value (ay) for each 
spin oy by calculating the product of the N 2 x 2 transfer matrices. Since some of the random 
matrix elements can be very small, floating point accuracy gives a lower limit of T = 0.05. 

First we address the scale-lengths of the equilibrium magnetization by computing the aver- 
age length [Z m ]av that separates two zeroes of the magnetization m(x). Figure ^ demonstrates 
how this length-scale changes with temperature, if we first scale away the T — 0-dependence 
on the field. A further collapse with the right combination of h r and T makes it possible to 
observe an universal scaling function for [Z m ] av 

Mav = [l d Uf(T/h 2 r / 3 ), (7) 

where the scaling function / — > 1 with T — > 0. The dependence of [Z m ] av on the combination 
of temperature and field strength does not follow an Imry-Ma-like scaling but is a consequence 
of entropic effects. The length [/ TO ] a v at a finite temperature is determined by both a zero- 
temperature scale ([i m ]av) and thermal fluctuations. The following argument can explain the 
scaling variable h 2 -^ 3 /T, analogously to spin glass chains in an external field Q. Once again 
consider the non-absorbing random walks which the domains consist of. Some of these inside 
a typical domain are such that the random walk sum of fields over the excursion is close to 
2 J. These almost-absorbing walks are the sequences (of spins) most likely to be flipped at 
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Fig. 5 - Mean Zeeman energy Hz — 2J corresponding to a particular domain length Id in a log-log plot. 
h r = 3.0, 2.0, 1.2, 1.0, 0.8, 0.6, 0.4, 0.2, 0.13, 0.08, 0.05 from top to bottom. The slopes of the straight 
lines are all within the interval 1.00 ± 0.05. The data are averaged over 10 disorder configurations. 
The straight lines represent least square fits. 

Fig. 6 - Scaling-plot for the average domain-lengths at finite temperatures: [l(T = 0)] a v/[{^T)]T versus 
T/£(h) for different values of the field strength h r , where the length scale oc hr . 



finite temperatures. The cost of flipping such a part of a domain is proportional to J, which if 
measured in terms of h r can be written with the help of the length-scale I of the non-absorbing 
excursions, h r ~ l/l. This is almost equal to the Zeeman-energy optimized over the excursion, 
which scales as Ef ~ h r l x l 2 . Equating the cost with the gain and solving for the energy scale 
(Ef) as a function of field gives rise to the Arrhenius factor Ef/T ~ hf^ 3 /T. 

As Fig. 1 demonstrates, the magnetization profile at T > differs from the GS due 
to domain wall fluctuations and internal cluster reversals. To study this quantitatively we 
introduce a parameter c S (0, 2) and define a reversal to be a sequence of spins for which 
\(<Ji)(T) — <Ji(T = 0)| > c holds; moreover the definition can be applied to both processes 
separately revealing interesting details. Since bulk reversal is always coupled with the breaking 
of two extra bonds one expects that domain wall fluctuations dominate. However, the former 
contributes a considerable portion to the total melting even at low temperatures (Fig. 0). The 
relative portion of bulk reversals at first grows with temperature for all values of c since the 
gain in entropy allows for more broken bonds. Moving the threshold c away from 1 and — 1 
respectively, a greater number of bulk segments are identified. Eventually for very large c 
even more bulk than boundary reversals are observed. The characteristic reversal rates are 
different for the two processes, and related via the empirical formula 

.Am. i^ m \ (o\ 

AT ' bulk = AT ' bound ' ^ ' 

with a w 1.63. Thus the change in magnetization with increasing T is stronger inside the GS 
domains than at their boundaries. These results are independent of the field strength. 

In conclusion, we have studied the magnetization properties of one-dimensional RFIM 
chains. These can be explained using of random walk arguments. While the GS structure 
is found to be a sequence of absorbing and non-absorbing excursions, the finite-temperature 
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Fig. 7 - Left: Fraction of melted spins as a function of temperature for different values of c. 
Right: Portion of those spins is displayed which reside inside bulk segments. L — 400, h r = 1.0. 

magnetization is complicated by thermal excitations. These are explained with the help 
of almost-non-absorbing walks. The results illustrate how a global optimization problem 
influences physics at T > in systems where the geometric arrangement of domain walls is 
crucial. Yet, extending the results to higher dimensions seems insolvable. 
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